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Abstract 

The object of this paper is a Bayesian analysis of the 
autoregressive model X t = PiX t _ t + £ t ,t=l,..,m and X t = 
P 2 -Xt-i + £ t ,t=m+l,..,n where 0 < ft, p 2 < 1, and e t is 
independent random variable with an exponential 
distribution with mean 9 1 but later it was found that there 
was a change in the process at some point of time m which is 
reflected in the sequence after s m is changed in mean 9 2 . The 
issue this study focused on is at what time and which point 
the change begins to occur. The estimators of m, ft, ft and 
9 1,9 2 are derived from Asymmetric loss functions namely 
Linex loss & General Entropy loss functions. Both the non- 
informative and informative priors are considered. The 
effects of prior consideration on Bayes estimates of change 
point are also studied. 
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Introduction 

In applications, time series data, the common issue, 
with characteristics which fail to be compatible with 
the usual assumption of linearity and Gaussian errors, 
for a variety of reasons. Bell and Smith (1986) 
concerned with the autoregressive process AR (1) X t = 
pXi_! + £j where 0 < [i < 1 and £/s are i.i.d. and non- 
negative, who considered the estimation and testing 
problem for three parametric models: Gaussian, 
uniform and exponential. For large series, non 
normality may not be of importance due to the 
additive nature of filtering processes; but otherwise for 
small series. Hence, model other than Gaussian, in 
particular, the exponential, is studied. This AR (1) 
model can be used as a model for water quality 
analysis. Let the initial level of pollutant beX 0 ; and 
random quantities be e lt e 2 , ■■■ Of that pollutants are 
"dumped" at regular fixed intervals into the relevant 
body of water and successive "dumping" a proportion 
(1 - p ) of the pollutant ( 0 < /? < 1) is "washed away". 
If ft- is level of pollutant at time t then ft= ftft^ + £,- 


where, £,- are assumed to have an exponential 
distribution. Some of the references for Bayesian 
estimation of parameters of AR (1) as defined above 
are A. Turkmann, M. A. (1990) and M. Ibazizen, H. 
Fellage (2003). 

Apart from AR (1) with exponential errors, the 
phenomenon of change point is also observed in 
several situations in water quality analysis. It may 
happen that at some point of time instability in the 
sequence of pollutant level of a river is observed. The 
issue this study focused on is at what time and which 
point the change begins to occur, which is called 
change point inference problem. Bayesian ideas, 
playing an important role in study of such change 
point problem has been often proposed as a valid 
alternative to classical estimation procedure 

A sequence of random variables ft, ...X m ,X m+1 ...X„ is 
said to have a change point at m (1 < m < n) if ft ~ 
ft ( x|ft ) (i=l,2...m) and Xi~ft(x|ft) (i=m+l,2...n), 
where ft ( x\ ft ft ft ( x\9 2 ). The situation is in 
consideration in which ft and ft have exponential 
form, but the change point, m, is unknown. The 
problem has also been discussed within a Bayesian 
framework by Chernoff & Zacks (1964), Kander & 
Zacks (1966), A.F.M. Smith (1975), P.N. Jani & Mayuri 
Pandya (1999), Pandya, M. and Jani, P.N. (2006), 
Pandya, M. and Jadav, P. (2009) and Ebrahimi and 
Ghosh (2001). The Monograph of Broemeling and 
Tsurumi (1987) on structural changes and a survey by 
Zacks (1983) are also useful references. 

In this paper, an AR (1) model with one change point 
has been proposed, where the error distribution is 
supposed to be the changing exponential distribution. 
In section 2, a change point model related to AR (1) 
with exponential is developed. In section 3, posterior 
densities of ft, ft, ft, ft ar| d m for this model are 
obtained. Bayes estimators of ft, ft , ft, ft and m are 
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derived from asymmetric loss functions and 
symmetric loss functions in section 4. A numerical 
study to illustrate the above technique on generated 
observations is presented in section 5. In section 6 , the 
sensitivity of the Bayes estimators of m when prior 
specifications deviate from the true values is studied. 

In this study, section 7 is about the simulation study 
implemented by the generation of 10,000 different 
random samples. The ending of the paper is about the 
detailed conclusion made on the study. 

ftX ) _ x + £ i; t = 1,2,3 
fi 2 x i-i + £ i . i — m + 

With X 0 fixed constant and, 0 < ft , ft < 1 . 


Proposed AR (1) Model 

Let { e n , 0 < s n , n > 1} be a random sequence having 
exponential distribution viz. 

The p. d. f. of the distribution is given by, 

/(£;lft) = l/6 1 .e~ £ ^ ei , i = 1,2,3, ... .m 

= 1 / 6*2 , i — m+1, n 

Further, let {Xn, Xn, >0, n > 1 } be a sequence of random 
variables defined as, 

. . m 

(2-D 

both cases, estimation will be same. 


In literature, there are two models for x 0 

Model A: x 0 is constant, x 0 =0 in particular; 

Model B: x 0 has the same distribution as that of e t 

Here, model A for its flexibility is examined, moreover, 
the likelihood function (conditional to x 0 ) for model B 
is exactly of the same form as that for model A. In the 

f(X 0 ,X 1 x n )= (±-) m .et-Sm + Pi s' m ) 

ft = YlU ft ; 

Since £; > 0, i = 1, 2, . . ., n. We have 

ft >ft.ft_ l7 and X t >0, i = 1, 2, ..., m, X t >ft. X;_ 1; X t >0, 
i = m+1, ..., n. and as a result 0 < /?i, /S 2 < 1 . 

Then, the likelihood function giving the sample 
information is obtained 

X=(X 1 ,X 2 X m ,X m+1 , ...X n ) is. 


(2.1) is the first order auto regressive process, AR (1), 
and m is unknown change point in the sequence to be 
estimated. For ft = ft = ft m = n . the model (2.1) 
reduces to the model studied by Bell and Smith (1986). 

Since, the AR (1) defined in (2.1) is Markov process, 
the joint p. d. f. of X Ql X l , X m ,X n is given by 

h ^ g/hbn- S m )/ s 2 _ e -(. s n~ s m)/ s 2 

c* V 

— Li = 1 * i - 1 

^(ft.ft. ft. ft, m | X_), with a prior joint density, say, 
^(ft.ft. ft. ft , m | X) representing uncertainty on the 
parameters values. 

^(ft.ft. ft. ft. 771 1 K) = 

\X) .g(P 1 ,p 2 A.&2.™ I jQ 

f fil Sp 2 W fff 2 LtPl.foAA.m I x) ■gtfl.PzAA.m \r)dp 1 dp 2 de 1 d& 2 


L(ft.ft, ft, ft.m I X) = . e~ A l e K (e 2 )-(n-m) e -n/e 2 

ft 

9i, 0 2 > 0, Xq — 0 

Where, 

A = s m -fts; 

B =S n -Sm-P 2 (Sn -S* m ) (2.2) 


Bayes Estimation 

The ML methods as well as other classical approaches 
are based only on the empirical information provided 
by the data. However, when there is some technical 
knowledge on the parameters of the distribution 
available, a Bayes procedure seems to an attractive 
inferential method. The Bayes procedure is based on a 
posterior density, say, 5 (ft, ft, ft, ft,m | X), which is 
proportional to the product of the likelihood function 


Using Inf ormative Priors for ft and ft, 9i and 61 

The first step in a Bayes analysis is the choice of the 
prior density on the distribution parameters. When 
technical information about the mechanism of process 
under consideration is accessible, then this information 
should be converted into a degree of belief on the 
distribution parameters (or function thereof). This 
problem is greatly aggravated when some 'physical' 
meaning can be attached to these parameters. 

It is also supposed that some information on these 
parameters is available, and that this technical 
knowledge can be given in terms of prior mean values 
pi, p 2 . Independent beta priors on ft and /i 2 with 
respective means /jl, fj2 and common standard 
deviation a viz are supposed. 

ftft)=[ft“ 1 - 1 (l- ft)* 1 " 1 ] I Pfti.ft), 0 

< 1 , Q.\ > 0 , > 0 
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9W 2) = [/?2 2-1 (l- ft)* 2-1 ] 1 P(a 2 ,&2), 

< /?2 ^ 1 / ^2 ^ ^2 ^ 0 


0 and#!, 6^. 


g(jn ) = 


If the prior information is given in terms of prior n — 1 

means pi, |_u and a, then the hyper parameters can be For simplicity, it is assumed that, a priori. 


obtained by solving 

a; = or 1 [C 1 - /h)/L 2 - M,cr £ ] i = 1,2. 

b t = p? (1 - i = 1,2. 

Let the joint prior density of 9 1 and 0 2 be 
g(^i' $) 


Pi,P 2 , @i> @i and m are independently distributed. The 
joint prior density of Pi,P 2 , 0\> @2 and rn say, 

9i(Pi, Pi, 0\, 0i,m) is. 


gi(Pi,Pi, ei.fy.m) =k 1 - 

Where, k 1 = 




P ( fl hl ) P ( <12*2 )(«-!) 


As in Broemeling et al.(1987), the marginal prior Joint posterior density of P\,P 2 , 0i> 02 and m say, 

distribution of m is supposed to be discrete uniform 9 \(P\,Pi, 0i> 0 2 ,m\ x) is, 

over the set { 1, 2, n - 1} and ndependent of pi, p2 

g 1 (Pt,p 2 , 0 y , 0 2 ,m | X) =ki% ( m+1 V^ e -B/b 

Pi 1 - 1 a - Pi? 1 - 1 p 2 2 -\i- p 2 ) h2 - x / h x {x) , 

hi (x) is the marginal density of X given as, where, 

hi(x) = k 1 Z n m -=\ him) , (3.2) 

Ti(m)= = rm r (tl — m). ((— )(5 m ) _m Appel iFi[ai,m, — b x , 1+ a 1> — ,1] ] 

l a l Sm J 

{(-h [(S n - Appel iFi[a 2 ,,n-m,-b 2 ,l + a 2 , S f^,l-\} , (3.3) 

1 a 2 S n -S m J 


Note-1 


Appel iFi(a,bi,b 2 ,c,x,y) is Appel Hypergeometric 
function defined as 


n r 

Appel 1 F 1 (a,b 1 ,b 2 ,c,x, y) = — — -J u a_1 (l — u) c_a_1 (l — ux)~ bl (l — uy')~ b2 du, 


Realja] > 0, Real (c-a) > 0. 

Integrating g x {Pi,P 2 , &i, 6^, m | X) with (P 1 ,P 2 ) and 
( 6\ , 0 2 ), leads to the posterior distribution of change 
point m. 

Marginal posterior density of change point m is given 


(3 - 4) 

T £ (m) is as given in (3.3). 

Marginal posterior densities of Pi and p 2 say, g t ( Pi I X) 
and g 2 {Pi IX) are as. 


9l(Pl I X) = ST=- 1 tie So 9 1 (Pv Pi , 0i, (h,m \ X) .d0 1 d9 2 dB 


=kiH n m -=\ (i - Pi) 


Ks» - 5 m )]- ( — } 


.Appel iFi [a 2 ,,n-m,-b 2l l + a 2 ? n _ S ' n , 1]| / ^(x) 

a n J 

gi (Pi i jd =fcizr=\ p?~k i - Pi ? 2 - 1 rmr J-~ m) {(£) (s m r m 

Appel iFi[ai,m,- bi , 1 + — ,1 ] ] / hi ( X ) 

s m ' 

Marginal posterior densities of and 0 2 , say cj\ ( 0] IX) and g t (0 2 IX) are given as, 

gi(0i\X) = xr=i Jo 1 Jo 1 f~gi (Pi, Pi, 01, 02, m I X)d6 2 d[3 1 d[3 2 
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=hr m -=\ ^ (m+1) d ^F(n ~ m){ (£) [(5 n - 

dppeZ 1 F 1 [a 2 ,,n-m,-ft,l + a 2 , ^~'! m ,l]/ ft© 


(3.7) 


And 


5i (ft I x ) = zr=\ Jo 1 Jo 1 Jr <71 (ft, ft, ft, ft,™ I X)d0 1 dp i d[3 2 

= ftzrft dfe {(ft) (s m )- 

Appel 1 F 1 [ a 1 ,m, — ft , 1+ a v — , 1] I / ft (v) 

5m J 

Using Non Informative Priors for |3i, ft, 0i and 02 variables, then the non-informative prior is. 

The non-informative prior is a density which adds no 92 (ft, ft, ft, ft,™ ) ix ^ 
information to that contained in the empirical data. If 


(3.8) 


Joint posterior density of 0i, 02, [3i, [32 and m say, 
92 (ft, ft, ft, ft,™ I ft) is given by 


there is no available information on ft, ft, ft, ft, m 
which are assumed to be a priori independent random 

= -L0-O+ 1 ) g-S m +/?iS£M grCn-m+l) e -S„+S m +/? 2 (S^-S^)/(% /f 2 (X) 

Where, ft (x) is the marginal density of x given by, 

— ft(m) (3.9) 

Where, 


T2(m)= FrnF(n — m) — : 


f (V - r (n “ m 5 + (V - - Sn 


(s ;xm) 


■(- 


(n-m)(s;-Sm) 


)] 


Marginal posterior density of change point m is given by 

n — 1 00 00 1 1 

(”'*>- I/JJJ 52 (ft, ft, ft, ft,™ I ft) d0 1 d0 2 d(3 1 dp 2 

(3.11) 


52 1 


m=l 0 0 0 0 
= T 2 (m) / x;; 1 , F 2 (m) 

Marginal posterior density of ft and ft are given as. 


5 2 (ft i K) = zrA flr (m+1) r ps;,M - 


. ( e 5 ™/® 1 — l) 


hh, - )~ (n ~ m > + (h, - - s- n +s - m )-( " ■ ~”0 

k (n-m)(s;-Sm) 


)/ ft© 


5 2 (ft I ft) 


fn — 1 ru _ ( n_m+1 ) e 


n-1 


(^n-O ■ 

e % 




Sm m + (Sm-Sm)~ 
(S^)(m) 


L ]/ ft© 


Marginal posterior density of ft and ft are given as, 

5^il|)=^aftrmF(n - m) _L_ 

ft (Sn - S m )-(" ~ m ft (5 n - S m - S’ +S^ )-fr ■ ~ m > M , s 

lv (n—m)(S2 —SL ) JJ/ftUJ 


5 2 (ft 1 I) = ftT ^=i r mr<JL ~ 


(n-m)(s;-s^) 

1 


-(n-m) 


[(- 


(•^n /?2 (^n *^m )) 

Sm m + (S m -S^r m 


( n— m) 


(S^)(m) 


-)] /ft© 


(3.10) 


(3.12) 


(3.13) 


(3.14) 


(3.15) 


Remark 1: For n=m, ft = ft, ft = ft the equations (3.5), densities of /? and 0 of AR(1) process without change 
(3.7) and (3.12), (3.14) reduce to the marginal posterior point. 
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Bayes Estimates under Symmetric and 
Asymmetric Loss Functions 

The Bayes estimate of a generic parameter (or function 
thereof) a based on a Squared Error Loss (SEL) 
function Li (a, d)= (a-d) 2 , where d is decision rule to 
estimate a which is the posterior mean. As a 
consequence, the SEL function is relative to an integer 
parameter, 

Li'(m, v) °° (m-v) 2 , m, v=0,l,2, 

Hence, the Bayesian estimate of an integer-valued 
parameter under the SEL function Li'(m, v) is no 
longer the posterior mean and can be obtained by 
numerically minimizing the corresponding posterior 
loss. Generally, such a Bayesian estimate is equal to the 
nearest integer value to the posterior mean. So, the 
nearest value to the posterior mean is regarded as 
Bayes Estimate. The Bayes estimators of m under SEL 
are the nearest integer values to (4.1) and (4.2), 


rn = YW=\ m T x (m) / ZV~=\ V(m) 

(4.1) 

m** = 

YW=i mT 2 (m ) / YW=\ T 2 (m) 

(4.2) 

Other Bayes estimators of a based on the loss functions 

L 2 (a, d) = 1 a - d 1 


t / (0, if | a — d | < e, e > 0 

L 3 (a, d) = j „ ' . 

( 1 , otherwise 


is the posterior median and 
respectively. 

posterior mode. 

Asymmetric Loss Function 


The Loss function L (a, d) provides 
financial consequences arising from 

a measure of the 
a wrong decision 


1 



rule d to estimate an unknown quantity a. The choice 
of the appropriate loss function is just dependent on 
financial considerations, but independent upon the 
used estimation procedure. In this section, Bayes 
estimator of change point m is derived from different 
asymmetric loss function using both prior 
considerations explained in section 3.1 and 3.2. A 
useful asymmetric loss, known as the Linex loss 
function was introduced by Varian (1975). Under the 
assumption that the minimal loss at d, the Linex loss 
function can be expressed as, 

Li (a, d) = exp. [qi (d- a)] - qi (d - a) - 1, qi # 0. 

The sign of the shape parameter qi reflects the 
deviation of the asymmetry, qi > 0 if over estimation is 
more serious than under estimation, and vice-versa, 
and the magnitude of qi reflects the degree of 
asymmetry. 

Minimizing expected loss function Em [Li (m, d)] and 
using posterior distribution (3.4) and (3.11), the Bayes 
estimate of m using Linex loss function is obtained by 
means of the nearest integer value to (4.3) under 
informative and non-informative prior , say, say m* L 

mV ■ 

m L — ~ y qi InHZ-J 1 e~^T^m)/r m -=\T^m)l 

mV = - 1 /q x In Em=\ T 2 T 2 (m)], (4.3) 

Where Ti (m) and T 2 {m) are as given in (3.3) and 
(3.10). 

Minimizing expected loss function Em [Li (|3i, d)[ and 
using posterior distributions (3.5) as well as (3.6), the 
Bayes estimators of pi and p 2 are obtained using Linex 
loss function as, 

[E e~ q i^i] 


“in [ £r=i Jo 1 d ft rmr ^~ vP [&, - 

Appel 1 F 1 [a 2 ,,n — m, — b 2 , 1 + a 2 , Sn — , 1 ]) / hi (10 

Sn — Sm J 

Pk == -- In [ > I 1 

Rl Jn 


1 - Vf 1 ^ 2 ” 1 ( 1 - y ? 2 ) i ’ 2-1 


(4.4) 


cL(3 2 rmfin — m) { (5 m ) 


771 = 1 


^ 71—771 + 1 

Appel 1 F 1 [ a.i,m, - b x , 1+ ,1] } / fhGD 


(4.5) 


s* 

Where Appel 1 F 1 [a 1 ,m, — b x , 1+ a Xl — 
Appel a 2 ,,n-m,- b 2 ,l + a 2 , n ~ 1 71 , 1] 
as in note 1 and fii(V) is as in (3.2). 


,1 ] | and 
are same 


Another loss function, called General Entropy loss 
function (GEL), proposed by Calabria and Pulcini 
(1996) is given by. 
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Ls(a, d) = (d/a)'' 3 - q 3 l n (d / a) - 1, 

Whose minimum occurs at d = a, minimizing 


using General Entropy loss function by means of the 
nearest integer value to (4.6) and (4.7) under 
informative and non-informative prior, say m* E and 

m* E * 


expectation E [Ls (m, d)] and using posterior density « 

m E 

gi(m | x ), i = 1,2. The Bayes estimate of in is achieved 

m* E = [£' 1 [m“ ( ' 3 ]] 1/q3 = m~ q3 T 3 (m)]“ 1/( ' 3 . 


m* E * = [E.im-^]] = E"=i m~^T 2 (m)/JZ "it T 2 (m)]- 1/<?3 , (4.7) 

Putting q3 = - 1 in (4.6) and (4.7), Bayes estimate, using posterior distributions (3.5) and (3.6), we get the 
posterior mean of m is acquired. Bayes estimates of (3i using General Entropy loss 

Minimizing expected loss function E [Ls (|3i, d)] and function respectively as, 

p; R = META rm r(n - m )^ ( S m r m ( ,,_^ )ll2 ^ Ks n - 


S * 1 

Appel 1 F 1 a 1 -q 3 m,-b 1 , l+a 1 -q 3 ,— ,1 

Appel 1 F 1 [a 2 ,,n-m,- b 2 ,l + a 2 , S f^Y-,l]} / h^X)} Uq3 (4.8) 

Jfi J 

p2E=kAr m ~=\ rm r(n -rn)-±- [(5 n - f {(S m )~ m 

a 2~ c l3 ai 

Appel i G [ a 2 - q 3 ,,n -m,- b 2 ,l + a 2 - q 3 , " _ m , 1 ]} 

Appel 1 F 1 | a!,!?!, — fq , 1 + .ljj/hi^)} 1/93 (4.9) 

Where Appel b 3 , 1+ a v -f- ,1]} and using posterior distributions (3.14) and (3.15), we 

s*_s^ S ’ n get the Bayes estimate P* E of |3i, i = 1, 2 using non 

Appel iF, a 2 ,n — m, — b 2> , 1 + a 2 , ,1 are same . r ,. . , ^ c . T 

1 11 z z s n -s m J informative priors, under General Entropy Loss 


as in note 1 and is as in (3.2). 


function is given as. 


Minimizing expected loss function E [Ls (|3i, d)[ and 

n** ft vn-1 i ■ [ f (5 n — 5 m ) m )+ (s n — s m — S n -t-S m ) 


p;; - fezsr.vm [(- 


(n-m)(s;-Sm) 


(■^m) 2 F 1 [ 1 — q 3 ,m, 2 — q 3 ) 1 /h 2 (x)| 

„ ( 1 S- m + (S m - s* )- m 

1 2E In — 1 /—i (Sm)(.rn) 


($n S m ) ^ m+1 ^ 2 F 3 [ 1 — q 3 ,n — m, 2 — q 3 , m ](1 — q 3 ) 1 / 1^2 (2i)l 


Where h 2 {X) is same as in (3.9). using posterior distributions (3.7) as well as (3. 12), we 

get the Bayes estimate 0 3E and (fp: of Qi using General 
Minimizing expected loss function E [Ls (0i, d)l and _ , T c 

° r L ' Entropy Loss function as. 


Z 1 1 

Hm + q 3 ) r(n - m) (5 m )^ m+ ^KS n - 5 m )] 

ai a 7 


-(n-m) 


Appel 1 F 1 \ a 1 ,m,-b 1 , 1+ a 1 ,| 2 -,l|] 4ppel 3 F 3 [ a 2 ,n - m, - b 2 ,, 1 + a 2 , Sn _ Sm ,l) / hQQ} 


Ce= {^ZViiFCn- m) 


r(-l+m+q 3 ) [ 


(n-m)(S*-S^) 


( 1 \-t+ m +93 / 1 \-l+m+q 3 1 I 

[-£) + t^) 


(4.13) 


s m - s m > 0,S m > 0,771 + q 3 > 0 , 
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Minimizing expected loss function E [Ls (02, d)] and the Bayes estimate 9 2E and &£ E of 02 using General 
using posterior distributions (3.8) and (3. 13), we get Entropy Loss function as, 

(?2 e = {fct 2T=- 1 rmr(n - m + q 3 )(±- ) [(5„ - (S m 

a 2 Vai/ 


Appel ^ a 2 ,n-m + q 3 ,~ b 2 ,, 1+ a 2 , " _ , l| 

dppel iFi [ «!, m, - h , 1 + a lt ,1 ] ] /h t Q0} 

( 1 V n n-l + rc-m + qs) 

^ = i^Z rm — — 


K- 4 r) 


— 1 +n— m+^3 


s- m + (s; - s^) 

(Sm)(m) 


( 

- s m - (5* - 5^) 

_ c* h— m x _ 1/ 3 

irH/*®} 


— 1+n— m+^3 1 


S m - Sm > 0 , s m > 0 ,m+ q 3 > 0 , 


Remark 2: For n=m, /?, = /? 2 , h, = ff 2 the equations (4.4), 
(4.8), (4.10) and (4.12), (4.13) reduce the Bayes 
estimates under asymmetric loss function of [3 and 6 of 
AR(1) process without change point under asymmetric 
loss functions. 

Numerical Study 

Let us consider AR (1) model as 

_ f 0.1 Xi_ t + , i = 1,2, ...12 

Xl - (0.8 X l _ 1 + £,- , i = 13,14, . . .20. 

Where, £ ( 's are independently distributed exponential 
distributions given in (2.1) with 

9i=0.6, 02=1. 20 random observations have been 
generated from proposed AR (1) model given in (2.2). 
The first twelve observations are from exponential 
distribution with 0i=O.6 and next eight are from 
exponential distribution with 02=1. pi and (3 2 
themselves are random observations from beta 
distributions with prior means pi = 0.1, 


\X 2 = 0.8 and common standard deviation o = 0.1 
respectively, resulting in ai = 0.001, 

b =0.09, a2=0.48 and b 2 = 0.12. These observations are 
given in Table-1. 

Posterior mean of m, 91, 91, (il , and (12 and the 
posterior median of m have been calculated.. Posterior 
mode appears to be a bad estimator of m. For a 
comparative purpose point of view, estimators under 
the non informative prior are also calculated. The 
results are shown in Table-2. The Bayes estimates 
m* L ,m* E of m, 9 1E , 9 2E of 9\ and 0 2 , P( E , P 2E of (11 and 
(I 2 respectively for the data given in Table-1 have been 
computed. As well as the Bayes estimates under the 
non informative prior and asymmetric loss functions 
which are also calculated, from which the result shown 
in Table-3 reveal that rn'i, m*E, P( E , p 2E , 0\ E and ff 2E 
are robust with respect to the change in the shape 
parameter of GE loss function. 


TABLE 1 GENERATED OBSERVATIONS FROM PROPOSED AR (1) MODEL. 
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TABLE 2 THE VALUES OF BAYES ESTIMATES OF CHANGE POINT 


Prior Density 

Bayes estimates of change point 

Bayes estimates of 
auto correlation 
coefficient. 

Bayes estimates of 
Exponential parameters 

Posterior Median 

Posterior Mean 

P> 


01 

02 

Informative 

12.02 

12 

0.1 

0.8 

0.6 

0.95 

Non informative 

12.22 

11 

0.12 

0.82 

0.63 

1.03 


TABLE 3 THE BAYES ESTIMATES USING ASYMMETRIC LOSS FUNCTIONS. 


Prior Density 

Shape 

parameter 

Bayes estimates of change 
point 

Bayes estimates under General Entropy 
Loss 


4 1 

q 3 

ml 

ml 

Pi e 

Pie 

ff\E 

&iE 

Informative 

0.8 

0.8 

12 

12 

0.83 

0.13 

0.13 

0.95 

1.2 

1.2 

12 

12 

0.82 

0.12 

0.12 

0.94 

1.5 

1.5 

12 

12 

0.81 

0.11 

0.11 

0.91 

Non 

informative 



m" 

m-E 

Pie 

PVe 

$1E 

$IE 

0.8 

0.8 

13 

13 

0.86 

0.16 

0.68 

1.4 

1.2 

1.2 

13 

13 

0.85 

0.15 

0.66 

1.3 

1.5 

1.5 

13 

13 

0.83 

0.13 

0.65 

1.2 


TABLE-4 POSTERIOR MEAN M* FOR THE DATA GIVEN IN TABLE-1 . 


F 1 

f.12 

m* 

m*E 

0.1 

0.6 

12 

12 

0.07 

0.8 

12 

12 

0.2 

0.4 

12 

12 


Sensitivity of Bayes Estimates 

In this section, the studied sensitivity of the Bayes 
estimates is obtained from section 4 with respect to 

change in the prior of the parameter. The means ]Ll \ , 

H2 and standard deviation o of beta prior have been 
used as prior information in computing the parameters 
ai, bi, a 2 , b 2 of the prior. Following Calabria and Pulcini 
(1996), it is also assumed that the prior information 
should be correct if the true value of Pi (P 2 ) is close to 
prior mean pi(|i 2 ) and is assumed to be wrong if Pi (P 2 ) 
is far from (j.i(pi 2 ). We have computed m* and m* E using 
(4.1) and (4.6) for the data given in Table-1 with 
common value of o =0.1. For q3 =0.9, considering 

different values of ( jLly , /1 2 ) ar| d result are shown in 
Table-4. 

The results shown in Table-4 lead to conclusion that m’ 
and m*E are robust with respect to the correct choice of 
the prior density of Pi (P 2 ) and a wrong choice of the 
prior density of Pi (P 2 ). Moreover, they are also robust 
with respect to the change in the shape parameter of 
GE loss function. 

Simulation Study 

In section 4, we have obtained Bayes estimates of m on 


the basis of the generated data given in Table-1 for 
given values of parameters. To justify the results, we 
have generated 10,000 different random samples with 
m=12, n=20, 0i=l.O, 02=2.0, (3i= 0.4, (32 = 0.6 and 
obtained the frequency distributions of posterior mean, 
median of m, iul’, m e* with the correct prior 
consideration. The results are shown in Table 5. We 
also obtain the frequency distributions of Bayes 
estimates of autoregressive coefficients given in section 
4 with the both prior considerations and the results are 
shown in Table 6 and table 7. The value of shape 
parameter of the general entropy loss and Linex loss is 
taken as 0.1. 

We also simulate several samples from AR (1) model 
explained section 2 with m=12, n=20, 0i =1.0, 0.6, 0.7; 02 
=2.0, 0.8, 0.9 and |3i=0.2, 0.4, 0.5; [3 2 =0.4, 0.6, 0.7. For 
each 0i, 02, (3i and [32, and Bayes estimators of change 
point m and autoregressive coefficients (3i and (32 using 
q3= 0.9 has been computed for different prior means pi 
and g 2 , Which lead to the same conclusion as for single 
sample, that the Bayes estimators posterior mean of m, 
and m*E are robust with respect to the correct choice of 
the prior specifications on Pi (P 2 ) and wrong choice of 
the prior specifications on P 2 (Pi) 
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Table 5: Frequency distributions of the Bayes estimates of the change point 


Bayes estimate 

% Frequency for j 


01-10 

11-13 

14-20 

Posterior mean 

18 

70 

12 

Posterior median 

10 

74 

16 

Posterior mode 

12 

73 

15 

ihl* 

13 

77 

10 

mE* 

14 

78 

8 


Table 6: frequency distributions of the bayes estimates of autoregressive coefficients b1 and b2 using general entropy loss function 

(informative) 


Bayes estimate 

% Frequency for 


0.1-0.3 

0.3-0.5 

0.5-0.7 

0.7-0.9 

Pie 

13 

82 

02 

03 

P 2 E 

11 

01 

84 

04 


Table 7: frequency distributions of the bayes estimates of (non-informative) 


Bayes estimate 

% Frequency for 


0.1-0.3 

0.3-0.5 

0.5-0.7 

0.7-0.9 

PlE 

13 

78 

06 

03 

Pie 

11 

01 

74 

14 


Conclusions 

Our numerical study shows that the Bayes estimators 
posterior mean of m, and mV. are robust with respect to 
the different choice of the prior specifications on (3i, (b 
and are sensitive in case prior specifications on both (3i, 
[32 deviate simultaneously from the true values. 
Numerical study also shows that posterior mean of m 
is sensitive when prior specifications on both Pi, P 2 
deviate simultaneously from the true values. 
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